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Abstract 

Laboratory experiments are conducted to examine granular wave patterns near on- 
set as a function of the container oscillation frequency / and amplitude A, layer 
depth -ff, and grain diameter D. The primary transition from a flat grain layer to 
standing waves occurs when the layer remains dilated after making contact with the 
container. With a flat layer and increasing dimensionless peak container acceleration 
r = An'^f'^A/g (g is the acceleration due to gravity), the wave transition occurs for 
r ~ 2.6, but with decreasing T the waves persist to F = 2.2. For 2.2 < F < 3.8, pat- 
terns are squares for / < fss and stripes for f > fss] H determines the square/stripe 
transition frequency fss = 0-33y^g/H. The dispersion relations for layers with vary- 
ing H collapse onto the curve X/H = 1.0 + l.l(/y^^)-i-32±o.03 ^Yien the peak 
container velocity v = li^Af exceeds a critical value, \gm ~ 3\/-D(7. Local collision 
pressure measurements suggest that v^m is associated with a transition in the hori- 
zontal grain mobility: for v > v^mj there is a hydrodynamic-like horizontal sloshing 
motion, while for v < v^^; the grains are essentially immobile and the stripe pattern 
apparently arises from a. bending of the granular layer. For y at ygm 

less than fss 

and V < Vgm, patterns are tenuous and disordered. 
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Fig. 1. Images of granular wave patterns near onset: (left) squares, / = 33 Hz 
(/* = 0.28), (center) squares and stripes, / = 37 Hz (/* = 0.31), and (right) 
stripes, / = 43 Hz (/* = 0.36). In each case I? = 0.17 mm bronze, T = 2.5, N = 4, 
container is SI, and the image size is 40.5 x 40.5 mm. 

1 Introduction 



Granular materials are collections of discrete solids for which even the simplest 
reahzations — ensembles of identical solid spheres interacting only via contact 
forces — exhibit a wealth of surprising behaviors [1]. Our work on granular 
media has focused on a dynamic phenomenon, the formation of subharmonic 
standing waves in vertically oscillated granular layers [2,3]. These strongly 
nonlinear waves form patterns of stripes, squares, hexagons and more complex 
patterns [4] , as well as localized structures called oscillons [5] , as a function of 
three dimensionless control parameters: the acceleration F, the frequency /* = 



f^H/g, and the layer thickness N — H/D. This paper concerns the square or 
stripe patterns that arise at the primary instability as F is increased; examples 
of these patterns are shown in Fig. 1. Despite being composed of discrete grains 
and having typical wavelengths A of only 20-30 D, the appearance of granular 
patterns is similar to that of patterns in fluid systems [6,7], the closest example 
being standing surface waves in a vertically oscillated liquid layer (the Faraday 
instability) [8]. 

Pattern formation in granular media can perhaps be described by continuum 
equations analogous to those used in fluid system [9-11], but the continuum 
equations remain largely untested by experiment. Aspects of granular patterns 
have been described using phenomenological models [12-14], but making the 
connection to real systems requires an understanding of granular media at a 
more microscopic level. Recently, molecular dynamics simulations have begun 
to yield new details concerning the behavior of granular waves [15]. In partic- 
ular, a simulation developed by Bizon et al. [16] quantitatively reproduces the 
granular patterns observed in three dimensions and thus allows calculation of 
quantities not readily accessible by experiment. 

In this paper, we will first describe our experimental system and then discuss 
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Fig. 2. (a) Experimental apparatus and (b) lighting. Patterns are illuminated from 
the side by light incident at low angles; high regions are bright and low regions are 
dark. 

some generic features of the flat-layer/wave transition. Next, we will present 
results for the scaling of the square/stripe transition frequency and the wave- 
length with particle size and layer depth. Last, we will use results from local 
measurements of the collision pressure and from dispersion relations to demon- 
strate and characterize the grain mobility transition. 



2 Experiment 

Our experimental apparatus is similar to that used in Refs. [3-5] and is de- 
scribed in detail in Ref. [17]. A layer of granular material is placed in the 
bottom of an evacuated container mounted on the armature of an electro- 
mechanical shaker and sinusoidally oscillated in the vertical direction (see 
Fig. 2(a)). Our granular layers consist of spherical non-cohesive and non- 
magnetic particles of bronze, 316 stainless steel, titanium, or lead with as- 
poured densities of 5.1, 2.6, 4.7 and 6.7 g/cm^ respectively, and with 0.05 < 
D < 0.8 mm and a poly-dispersity of ±10%. Data is obtained for D = 0.17 



3 



Table 1 

Container specifications. 

Container Shape Size (mm) Base (mm) Mass (g) 



S2 



SI 



C 



circular 126 (dia.) 
square 162 (side) 
square 106 (side) 



12.7 



12.7 



10.2 



2170 



1150 



1460 



mm bronze particles except as noted. Three different containers are used (see 
Table 1); each has a polished aluminum base to prevent static charge accumu- 
lation and clear plastic side walls and top for lighting and visualization respec- 
tively. An accelerometer mounted on the underside of the container measures 
the acceleration. The patterns are visualized using as a light source a ring of 
strobed light-emitting diodes, which encircle the container and illuminate the 
layer from the side (Fig. 2(b)). Images are acquired by a digital CCD camera 
mounted on axis. 

The momentum transfer generated during layer-container contact is measured 
by an acceleration compensated pressure sensor with a resonant frequency of 
300 kHz (PCB 112 A22). This sensor has a circular sensing area 5.5 mm in 
diameter; the sensor size can be compared with the wavelength of the pattern, 
which varies from 2 to 40 mm for the present measurements, and the particle 
size, typically 0.17 mm. The sensor is flush mounted (±0.02 mm) 13 mm from 
one container sidewall and 36 mm from the other in the bottom of cell S2. 
To avoid high sampling frequencies, a circuit consisting of a peak detector, a 
level shifter, and a sample-and-hold is used to record the maximum pressure 
from each collision of the layer with the container. 



3 Wave Onset 

Here we investigate the transition from a flat layer to standing waves in the 
regime where the layer free- flight time t fn is less than one period of oscillation 
of the container T. The layer state is characterized by local collision pressure 
measurements [18]. Related work is reported in Ref. [17] and in Rcfs. [19,20], 
where in the latter, global force and optical measurements characterize the flat 
layer state. For F < 1, the layer is always in contact with the plate and there 
is no significant relative grain motion. For F > 1, we identify three distinct 
stages of layer motion during each cycle: free-fiight - layer not in contact with 
container; impact - layer and container collide; contact - layer and container 
in contact and moving with the same vertical velocity. Impact imparts relative 
kinetic energy to the grains which, if large enough, enables the layer to change 
its configuration. 
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Fig. 3. Non-dimensionalized layer pressure for (a) T = 1.5 (flat layer) and (b) T = 3.1 
(waves) with /* = 0.48 (/ = 40 Hz) and = 9. Dashed curve is a plot of a{t)/g + l. 
The insets show expanded portions of the collision region and illustrate the increase 
in collision time associated with waves onset. 

We first describe some generic features of the layer-container pressure in the 
vicinity of the flat-layer/waves transition. A time series for the local pressure 
exerted by the granular layer on the container below the onset of the waves 
is plotted in Fig. 3(a). The pressure is expressed in non-dimensional form, 
P* = P/pgH, where P is the pressure and p is the as-poured layer density. 
P* = indicates the layer is in free-flight, the sharp jump in P* occurs at 
impact, while contact is shown by P* ~ o,{t)/g + 1 — the pressure calculated 
using the measured plate acceleration and assuming a solid layer attached to 
the plate (dashed curve). Several features of P* should be noted. As the inset 
in Fig. 3(a) emphasizes, the collision duration is short, less than 0.02T, which 
indicates the layer is relatively compact at the time of coUision. Furthermore, 
immediately following the collision, P* = a{t)/g + l, which shows the collision 
is strongly inelastic (if the layer bounced, P* would be reduced). Finally, P* 
and a(t)lg + 1 go to zero simultaneously, which indicates the entire layer is 
moving upward with the container velocity when free-flight begins. 

To contrast the layer behavior above and below wave onset, a time series of 
P* in the wave regime is presented in Fig. 3(b). Here the layer motion consists 
of just two stages, a free-flight followed by an extended coUision. There is no 
longer a contact stage as indicated by the observation that P* is always larger 
than a(t)/g + 1, and that a{t)/g + 1 goes to zero before P* goes to zero. From 
these observations, we infer that the layer is expanded at impact, that at no 
time during the collision is the entire layer totally at rest with respect to the 
container, and that the coUision continues for a short time after a{t) becomes 
less than g. 

At the onset of waves, the maximum collision pressure per cycle, P^ax^ 
creases rapidly, as Fig. 4(a) illustrates. This decrease is associated with an 
increase in layer dilation (see Fig. 3). With increasing F, the decrease in P^^x 
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Fig. 4. The hysteresis in the transition from a flat layer to waves is revealed in 
measurements of (a) maximum collision pressure and (b) dimensionless flight time 
versus acceleration amplitude for increasing (solid-line) and decreasing (dashed-line) 
r (/ = 26 Hz (/* = 0.31), N = 8.2). The inset in (b) plots the difference At 
between the flight time calculated from the inelastic layer model (o in (b)) and the 
experimental flight time. Experimental data are averaged into 0.01 intervals in F 
and averaged over approximately 20 full ramps of T. 

occurs at r = 2.7. However, the transition is hysteretic and the layer reverts 
to the flat state at F = 2.2, where P^ax increases abruptly (see Fig. 4(a)). 
For the flat layer, P^g^ increases with increasing F, while in the wave state, 
Pmax decreases with increasing F [21]. Additional measurements using steel 
and titanium particles show a somewhat different behavior [17]. Before the 
transition to waves (F 1.9), P* in these larger restitution coefficient mate- 
rials drops from 80 to approximately 40 but then continues to increase until 
the onset of waves at which point P* drops to 20 as for the bronze layers. 

Figure 4(b) compares the layer free-flight time measured from the experiment 
with the time calculated from the completely inelastic layer model [22,3]. When 
the layer is fiat, the measured and calculated flight times are nearly identical. 
However, when waves arise, the measured flight time is smaller than the model 
predicts (see inset in Fig. 4(b)). This discrepancy arises because the velocity of 
a portion of the layer is less than the container velocity when a{t) becomes less 
than g. These slower grains reduce the effective layer take-off velocity, which 
decreases tfu. For increasing F, there is a sudden decrease in tfu when waves 
arise. However, for decreasing F, and unlike P^axi ^fit continuously decreases 
until it equals the model value at which point the layer reverts to the fiat state. 
The equality of measured and model values of tfu at this transition implies 
that for waves to exist, a portion of the layer must be dilated when the bottom 
of the layer leaves the container. This necessary condition appears to be due 
to the strongly dissipative nature of the granular layer, which rapidly removes 
kinetic energy when layers remain in contact with the container. 

The qualitative features of the transition from a fiat layer to waves, shown in 
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Fig. 5. Stability diagram showing pattern transitions with increasing (o) and de- 
creasing (A) r as a function of (a) /* with dimensionlcss layer thickness = 9, 
and (b) A'^ with /* = 0.36 (/ = 30 Hz). Hysteresis decreases with increasing /* 
and increases with increasing N. In (b) for A?^ < 2 the hysteresis vanishes; as F is 
increased in this regime the layer dilates into a uniform low density gas-like state. 

Figs. 3 and 4, are the same for the entire range of H, D, and / we examine. 
However, other features of the transition are non-trivial functions of H and 
/. For instance, the F value at which the transition from a fiat layer to waves 
occurs, Fc, is plotted in Fig. 5(a) as a function of /*. For increasing F, Fc| 
decreases monotonically with /*, while for decreasing F, Fc^ is nearly constant 
except for a shght local increase near /* 0.4 [23]. The observation that the 
existence of waves depends on the layer never reaching the contact stage gives 
the following necessary condition for wave existence 

T — tfit — tcoll — h/Vcoll, (1) 



where h is the layer dilation (actual layer thickness minus H), tcoii is the 
collision duration, and Vcoii is the collision velocity. At the flat- layer /waves 
transition for decreasing F, F^i is nearly constant which implies h oc /^^ since 
T — tfit and Vcou arc both proportional to f~^. For the data in Fig. 5(a), 
Fc| ~ 2.2 which gives h = 0.25vcoiitfit- Similarly, applying Eq. (1) in the flat 
layer regime where Fci is a function of /, yields h oc f"^'^ — a slower decrease 
with increasing / than for waves. 

Fc is also a function of A^, as Fig. 5(b) shows. T^i is nearly independent of 
A^, while Fc| is roughly proportional to A^. The transition as measured by the 
rapid drop in Pmax is no longer subcritical ior N < 2 — the same A^" value 
below which coherent wave patterns disappear. Interpreting these results in 
light of Eq. (1), the increase in F^i with increasing A^ is due to the greater 
dissipation in the thicker layer which reduces h. The disappearance of waves 
for A^ < 2 is likely a result of large grain velocity fluctuations that destroy 
coherent layer motion. This interpretation is also supported by experiments 
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with different materials [17]: materials with a restitution coefficient larger than 
that of bronze, such as stainless steel and titanium, require deeper layers than 
does bronze to generate waves. Lead, which has a smaller restitution coefficient 
than bronze, requires a layer of only N — 1.1 for waves to exist. Additionally, 
in layers deep enough to form waves, we find that higher (lower) restitution 
materials have higher (lower) values of Fc for equal N. 

In Ref. [16], Bizon et al. note that waves in their simulations form when the 
time for peak growth is greater than that for peak decay. Prom this observa- 
tion, they calculate ~ 2.47 independent of either /, D, or H. The data 
presented above is in reasonable agreement with their prediction. Their obser- 
vation does not address the internal layer state and thus can not predict the 
value of and its dependence on / and N . With results from simulations, 
it will be possible to measure h and thus check the validity of the criteria for 
wave onset proposed in Eq. 1. 



4 Square/Stripe Transition 

The wave patterns at the primary bifurcation are squares for low /* and 
stripes for high /*, while for intermediate /*, both patterns arc observed si- 
multaneously (see Fig. 1). To quantify the square/stripe transition, we divide 
the pattern images into 3A x 3A regions, calculate the spatial power spectra, 
radially average the power within one full width at half maximum of the dom- 
inant wavenumber, and subtract the mean to obtain I{9). The autocorrelation 
function C{ti/2) = {I{6)I{6 + 7r/2))g/ {P{9))q is then used to characterize the 
pattern. For perfect squares C{ti/2) = 1, whereas for an image of perfect 
stripes C{tt/2) ~ 0. Differences between the measured and ideal values {i.e., 
C(7r/2) < 0) are due to the finite angular width of the spectral peak. Fig- 
ure 6 is a plot of C{tx/2) versus /*, which shows the transition from squares 
to stripes occurs in a narrow frequency range about the transition frequency 
f*g — 0.31, with C(7r/2) nearly constant above and below /*^. The value of f*g 
is not sensitive to details of the method: plotting I{9max =t^/2), where O^ax is 
the azimuthal location of the maximum value of I{0), as well as varying the 
sub- image size, yield the same result. 

Figure 7(a) is a plot of f*g as a function of F for fixed N . There is little 
change in /*^ as F is increased from near wave onset to close to the tran- 
sition to hexagons [3]. Figure 7(b) plots /*^ against N. Although N varies 
by more than an order of magnitude, the non-dimensionalized square/stripe 
transition frequency remains nearly fixed. Included in Fig. 7(b) are additional 
measurements that indicate f*g also does not depend strongly on D, material, 
or container size or shape. 
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Fig. 6. The transition from square patterns at low frequencies to stripe patterns 
at high frequencies is indicated by the decrease in the autocorrelation function of 
radially averaged power spectral density at 6* = 7r/2 (F = 2.5, A'' = 4, SI). The 
solid line is a fit to a + 6tanh[c(/* — f*g)] with f*g = 0.31. The labeled /* values 
correspond to the images in Fig. 1. 

Bizon et al. have proposed that the square/stripe transition occurs when the 
distance that the layer falls from its maximum height to the container is equal 
to the layer depth: for fall distances larger (smaller) than H, squares (stripes) 
are the preferred pattern [16]. Using a simple model that considers the layer 
as a totally inelastic object with no internal modes [22,3], they calculate f*^ ~ 
l/VS ^ 0.35. The data presented here gives an average value of f*^ ^ 0.33 
over a large range of experimental parameters. The preceding measurements 
as well as the criterion of Bizon et al. suggest that the square/stripe transition 



0.4 



0.3 



0.2 



0.1 



(a) 



0.4 



0.3 



0.2 



0.1 



(b) 



o 

• + 

□ ^ ^ o o o o 



2.5 



3.5 



10 



Fig. 7. (a) Square/stripe transition frequency as a function of F for N = 13 (con- 
tainer SI), (b) Square/stripe transition frequency as a function of N: O (T = 3.0, 
0.17 mm, C), + (F = 2.5, 0.17 mm, 5"^), x (F = 3.1, 0.66 mm, SI), □ (F = 2.5, 
0.46 mm lead, C), * (F = 3.0, 0.07 mm, C), • (F = 3.0, 0.78 mm, C), O (F = 2.8, 
0.17 mm steel, C), A (F = 2.8, 0.17 mm titanium, C). 
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is primarily a function of the layer depth H. Layers with differing and D but 
equal H make the transition from squares to stripes at the same / (except for 
layers with N ~ 0(1) where the pattern weakens due to large relative particle 
velocities). This result implies that the layer acts as a continuum, i.e., the 
size of the particles does not affect the bulk behavior although the discrete 
nature does. In contrast, the dependence of the critical acceleration amplitude 
for waves on and /, where scaling with H fails (Fig. 5(b)), suggests that 
the physics determining subcriticality and square/stripe pattern selection are 
unrelated. 



5 Wavelength Scaling 



The dispersion relation for surface waves in inviscid fluids in the absence of 
surface tension and in the shallow water limit {X > H) is X/ H — f'^^Jg/H. 

When plotted in this form, which uses H and '^H/g as the characteristic 
length and time scales respectively, the dispersion relations for different depth 
layers lie on a single curve. In granular layers, however, there are two natural 
length scales - layer depth and grain diameter. We show in this section that 
below a critical frequency, wavelength scaling is governed by H as is the case 
for fluids, while above this frequency, scaling with H fails. In the next section 
(Sec. 6) , we will show that the scaling breakdown is associated with a transition 
in the horizontal grain mobility. 

Figure 8(a) presents results for the wavelength as a function of frequency 
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Fig. 8. Wavelength as a function of frequency for layers of constant dimensionless 
depth N = H/D and varying particle diameter D. (a) A versus / for 0.07 mm (o), 
0.08 mm (□), 0.14 mm (O), 0.17 mm (A), 0.23 mm (<), 0.33 mm (v), 0.66 mm (I>), 
and 0.78 mm (+) (F = 3.0, N = 5, SI), (b) Same data plotted with dimensionless 
variables using H as the fundamental length scale. Near /* = 0.4 in (b), there is a 
kink in the dispersion curve, which moves to higher /* as N is increased. 
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Fig. 9. Dispersion relation for layers of constant particle diameter and varying di- 
mensionless depth, N = H/D: (a) Wavelength as a function of frequency: = 2.2 
(o), N = 3.1 (□), N = 4.0 (O), N = 5.4 (A), (<) N = 6.8, TV = 8.7 (v), N = 10.7 
(0), = 13.8 (+), AT = 18.1 (x), N = 22.8 (*), and AT = 30.7 (•) (D = 0.17 
mm, r = 3.0, container C). (b) Same data plotted with dimensionlcss variables. 
The dimensionlcss data in (b) are fit to Eq. 2; the inset shows the residuals. As N 
increases, the /* range over which data collapse occurs increases. 

for layers of bronze particles with N = 5, T = 3.0, and a range of particle 
sizes (0.07 < D < 0.8 mm). The wavelength decreases with increasing / 
and increases with increasing D. When H is used to non-dimensionalize the 
wavelength X* — X/H and the frequency /* = f-sju/g, the data collapse onto 
a single curve (Fig. 8(b)). A well-defined kink in the dispersion curve occurs 
near /* = 0.4, above which A* decreases more slowly with increasing /*. A 
similar slowing has been observed in other experiments but only for one [24] 
and two [16] disjoint values of D. 



Because is constant for the data presented in Fig. 8, it is unclear whether 
H or D is the length scale governing the collapse since H = ND. To test 
whether the data collapse obtained in Fig. 8(b) is due to scaling with H, we 
plot in Fig. 9 the wavelength versus frequency with fixed D — 0.17 mm but 
with A'^ varying from 2 to 31. A increases with increasing H and again the 
unsealed data collapse when A* is plotted as a function of /*. Similar results 
were obtained by Bizon et al. but only for two values oi N [N = 2.7, 5.4) [16]. 

The rescaled data in Fig. 9(b) fit the functional form X* — a -\- hf*'^ with the 
best fit given by 

A* = 1.0 + l.l/*"^-^^^°-°l (2) 



The residuals from the fit to the rescaled dispersion data in Fig. 9(b) (see inset) 
for /* > 0.4 exhibit a marked decrease in the quality of the data collapse, which 
is a function of H and which is characterized by a slowing in the decrease of 
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A* with increasing /*. This slowing is more pronounced in small N layers. We 
will return to this observation in the next section. 

Three points should be noted concerning the dispersion relation data. First, 
although the granular data collapse below the kink in the dispersion curves 
when scaled with H as for fluid surface waves, the dispersion relations for 
these two media are different. Second, the exponent in Eq. 2 is not equal to 
—2 as was found in our earlier work on glass beads in air [2], by Metcalf 
et al. [24] for glass beads in an evacuated container, and by Clement et al. 
in a two-dimensional layer of aluminum balls [25]. Possible explanations for 
the larger value of exponent found in the other studies include the influence 
of air viscosity, static charging in non-conducting grains, small aspect ratios, 
and drag associated with the side walls. Finally, the fit predicts a minimum 
wavelength equal to H for large /. The A* = 1 limit is clearly unattainable 
in small N layers due to the finite size of the layer's constituent grains; thus, 
scaling with H will fail for small N. However, as we will show next, breakdown 
of scaling normally occurs before this limit is reached due to a reduction in 
the grain mobility. 



6 Grain Mobility Transition 

As discussed in the previous two sections, the layer depth is the characteris- 
tic length scale which determines the scaling properties of the square/stripe 
transition and the wavelength. However, as mentioned in Sees. 1 and 5, the 
grain diameter D is also expected to influence the layer response. Consider 
the initial velocity required to raise a grain a distance D in the presence of 
gravity, vd = \f2gD. If the grain is embedded in a plane of identical particles, 
horizontal motion will only be possible iivo > \f^WD\ for lower velocities, the 
grain will be stuck. Assuming further that the fluctuation velocity of a typical 
grain with respect to its neighbors is proportional to the peak container ve- 
locity 27rAf, a relevant dimensionless parameter characterizing grain mobility 
is V = 2T{Af / \fDg. Section 6.1 presents evidence of a qualitative change in 
the layer response at a particular value of v which we call the grain mobility 
transition, v^^ ^ 2.5, for large F or small N. Section 6.2 then shows that this 
change appears related to the loss of horizontal grain mobility for v < v^^. 

6. 1 Changes in Layer Response 

The value of /* below which data collapse occurs for constant D and varying 
is marked by a kink in the dispersion curves (see Fig. 9(b)) and is an increasing 
function of A^. The data in Fig. 8(b) for constant A^ and varying D show a 
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Fig. 10. Dimensionless container velocity at which A* begins to decreases slower 
than Eq. 2, as a function of N for the data in Fig. 9(b) with F = 3.0 (o) and for 
additional data [17] with identical particles but with F = 2.5 (A). For F = 3.0, Vgm 
is nearly constant, while for F = 2.5, Vgm increases gradually with increasing N. 

similar kink which is, in contrast to the data for constant D, independent of 
H for fixed A^. To characterize the location of the kink. Fig. 10 presents a plot 
of V at the kink {e.g. Vgm) as a function of for F = 3.0 and F = 2.5 [26]. For 
r = 3.0, Ygm ~ 2.5 and shows no systematic dependence on A^. This is nearly 
the same as ygm — 2.6, the value found for varying Z), constant A^ — 5 and 
r — 3.0 (see Fig. 8). For F = 2.5, Vgm is slightly larger and is a slowly (slower 
than y/N) increasing function of N. 

To summarize, for v > Ygm, A* is independent of H and scales with /*, while 
for V < Ygm: A* scahug with H fails. This breakdown appears to be related 
to a transition in the layer response when v is small. Observations aX v ^ 0.5 
using relatively large grains show that there is little if any horizontal particle 
motion at the surface. In contrast, for v > Vg^, grains slosh back and forth - 
in thinner layers (A^ ~ 4) the sloshing is so vigorous that the layer depth goes 
to zero in the pattern minima. To illustrate the difference between these two 
wave regimes. Fig. 11 presents a plot of A* versus v for A^ = 5 and F = 2.5 
and F = 3.0. The wavelength is approximately 10 percent larger for F = 2.5 
than for F = 3.0 when v > Vg^ ~ 2.6, but for v < Vgm, A* for F = 2.5 jumps 
to 1.5 times A* for F = 3.0. If A* is instead plotted versus /*, the change in 
behavior is even more evident: for /* below the transition A* increases with 
increasing F, while for /* above. A* decreases with increasing F. 

In addition to the breakdown in A scaling with H, a qualitative change in 
the collision pressure is also observed for v ^ 3. The v dependence of the 
F averaged coUision pressure {P^ax)r = ItZ7 ^*(^)^^ plotted in Fig. 12 
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Fig. 11. A* = X/H versus v = 2tt Af/./D^ for T = 2.5 (•) and T = 3.0 (o) (0.46 mm 
lead, N = 5, container C). Inset: The A ratio jumps for v < 2.6, which indicates the 
grain mobiUty transition. 

for both the flat-layer and wave regimes. For the flat layer, {Pmax)r increases 
continuously with increasing v. However, in the wave regime {P^ax)r nearly 
constant for v > 3 but decreases in a similar fashion to the flat layer for v < 3 
[27]. The latter result suggests that the flat layer and wave states are similar 
for V < 3 as would be expected below the grain mobility transition. 

Finally, as Sec. 4 describes, the frequency of the square stripe transition is 
given by fgg — 0.33^g/H and shows no systematic dependence on the grain 

size. However, in thin layers when Ty/N < 2Ti\gmfssi the frequency associ- 
ated with the grain mobility transition fgm is less than the frequency of the 
square/stripe transition. In this case and for v < v^m, the resulting wave 
patterns are noticably more tenuous and have significantly shorter spatial 
correlation lengths (on the order of A) than do the corresponding patterns for 

fss ^ fgm- 



6.2 Loss of Horizontal Grain Motion 



In this section, the breakdown in A scaling with H and the decrease in the 
collision pressure arc identified with a transition in the horizontal grain mo- 
bility. We present data suggesting that for v > Vgm, the local layer height is 
changed by lateral grain motion, while for v < Vgm, grains are essentially im- 
mobile [28]. In Ref. [20], Mujica and Melo similarly propose that waves at low 
container velocities result from bending of the layer and not from horizontal 
grain transfer. Also, in our direct visual observations we observe significant 
horizontal particle motion for waves at large v but no noticeable horizontal 
grain motion for waves at small v. In general for waves at both high and 
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Fig. 12. F-averaged collision pressure versus container velocity below wave onset for 
1.2 < F < 1.8 (o) and above wave onset for 2.8 < F < 3.2 (•) {N = 9). For v > 3, 

{Pmax)r ^^'^ wave regime is constant, while for v < 3, {Pmax)r ^'^^ waves and 
flat-layer show similar behavior. The line at {Pmax)^ — ^ represents the calculated 
minimum value for P^^x at F = 3.0. 

low V, patterns persist when the container is rapidly brought to rest. When 
the stationary container is lightly and repeatedly tapped, the patterns at low 
V slowly disappear without any apparent grain motion on the layer surface, 
whereas the high v patterns exhibit significant grain rearrangement. 

To quantitatively investigate the proposed grain mobility transition, fluctua- 
tions in the local collision pressure are studied. In the presence of waves and 
for v > Ygm, relatively large fluctuations in the collision pressure associated 
with changes in the local layer height due to grain motion are expected. Con- 
versely, for V < Vgrn fluctuations should be smaller because the layer height is 
essentially constant. Figure 13(a) plots the relative standard deviation of the 



maximum pressure (Jp^^^ — J {l\PmaxY/ {Pmax) versus v for both the flat-layer 



and wave regimes where ^J (APmax)"^ and {Pmax) are the standard deviation 
and mean of Pmax respectively. For the flat-layer, crp^^^ is small and nearly 
independent of v. For the wave regime, crp^^^ shows three distinct behaviors — 
v < 3: (Tp^^^ is independent of v and equal to the flat-layer value; 3 < v < 7: 
^Pmax increases with increasing v and is larger than the flat-layer value; v > 7: 
^Pmax decreases with increasing v. 

As a further check that the equality of (Jp^^^ for the flat-layer and waves 
for V < 3 is due to a loss of horizontal grain motion. Fig. 13(b) presents 
measurements of the autocorrelation of Pmax at a delay of 4T, Cp^^^(4T) [29]. 
Cp^^^{4T) is sensitive to periodic variations in Pmax even when the intrinsic 
pressure noise is larger than the fluctuations associated with grain motion. As 
is true for (Jp^^^, Cp^^^^AT) for the flat-layer is independent of v. For waves, 
the dependence of Cp^^^{AT) on v can also be divided into three regimes — 




15 



0.8 



-'P 

max 

0.4 



(a) 



^ 0.5 



4 ~ 6 



10 



(b) 



a » 



4 ~ 6 



10 



Fig. 13. Grain mobility below wave onset at F = 1.9 (o) and in wave regime at 
F = 3.0 (•), characterized by (a) pressure fluctuations and (b) autocorrelation of 
peak pressure at a delay of 4T (iV = 9). Near v = 3, both quantities change, which 
marks the grain mobility transition. 



V < 3 : Cp^^^(4T) is constant and equal to the flat- layer value; 3 < v < 7 : 
Cp^^^{AT) is constant and larger than the flat-layer value; v > 7 : Cp^^^{AT) 
decreases with increasing v. The transition in Cp^^^(4T) to the flat- layer value 
is sharp and occurs at the same v where (Tp„„^ for waves reaches its minimum 
value. Also, note that (Jp^^^ peaks at the velocity (v = 7) at which Cp^^^{4T) 
begins to decrease for increasing v. 

Our interpretation of the three regimes for (7p„^^ and Cp^^^ shown in Fig. 13 
is as follows. For v < 3 horizontal grain mobility is strongly reduced and 
the layer depth is everywhere equal. Waves in this regime are likely bending 
waves [20] since mass transfer, dilational waves [19], or any other mode giving 
rise to a periodic variation in pressure would produce positive correlations 
in the pressure fluctuations. For v > 3 waves are mass transfer waves. We 
speculate that the decrease in (Tp^^^ and Cp^^^ above v = 7 results from 
the disordered patterns found in this regime which possibly result from an 
instability associated with the more rapid growth of the wave amplitude in 
comparison to A as / decreases [27]. 

Having made the case for a grain mobility transition in terms of crp„„^ and 
CprnaxJ '^^ ^o'^ examine the dependence of the transition on v, N, and D. 
Figure 14 presents crp^^^ for varying F and nine distinct / values ranging from 

as a function of v. Most importantly. Fig. 14 shows 
that variations in F at constant / produce the same changes in crp^^^ as do 
variations in / at constant F (see Fig. 13(b)). This finding strengthens our 

proposition that v = V / {2-k f ^ D / g) is the correct dimensionless parameter 
specifying the grain mobility. The data also show that v^^ ~ 3, which is 
in agreement with the value obtained from the pressure fluctuation data for 
F = 3.0. 
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Fig. 14. Normalized collision pressure fluctuation in wave regime versus dimension- 
less peak container velocity for varying r and f^/Dj^ = 0.05 (•), 0.07 (□) , 0.09 
(filled O), 0.11 (A), 0.14 (filled <), 0.17 (v), 0.21 (filled >), 0.25 (o), and 0.33 (filled 
□) (AT = 9). Note the similarity with Fig. 13(a). 

Figure 15(a) examines the dependence of grain mobility on layer depth by 
plotting Cp^^^{4T) versus for v = 3.4 and v = 4.2. By choosing v slightly 
above Ygm ~ 3 (for N — Cp^^^{AT) is sensitive to the grain mobility transi- 
tion. Below = 14, the layer response for the different v is similar. For N < 5, 
Cp^^^(4T) decreases to near at ^ 2, while for 5 < < 14, Cp^^^(4T) 
decreases slightly with increasing A^. The dependence of Cp^^^{4T) on A^ at 
small N appears similar to that of C'p^^^(4T) on v at large v, suggesting that 
the two limits of large v and small N are related. For N > 14, Cp^^^{AT) 
continues to decreases slowly with increasing A^ for v = 4.2, but, for v = 3.4 
and > 14, Cp^^^(4T) drops to approximately the flat layer value. Possibly, 
the reduction in grain mobility with increasing A^ is due to a decrease in grain 
velocity associated with an increase in the grain collision frequency. A more 
mundane explanation is as N increases, the layer mass becomes significant 
with respect to the container mass. This reduces the change in layer velocity 
at impact and subsequently decreases the effective driving. 

Finally, Fig. 15(b) looks at the effect of varying D on the grain mobihty by 
plotting Cp^^^(4T') versus v for D = 0.17, 0.33, 0.46, and 0.66 (pressure data 
for constant v and varying A^ was not collected). The dependence of Cp^^^^^{'iT) 
on V and the value of Ygm ~ 3 is in accord with the data for constant D and 
N and varying / shown in Fig. 13 and is consistent with our proposal that 
Ygm is independent of D. In addition to bronze. Fig. 15(b) also includes data 
for layers of lead and steel particles. 

In summary, Ygm values obtained from local pressure measurements and from 
identifying the breakdown in wavelength scaling with H are in good qualita- 
tive agreement. Both show the same weak dependence on A^ for small F and 
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Fig. 15. (a) Autocorrelation of peak pressure versus N showing grain mobility tran- 
sition at A'^ = 14 for v = 3.4 (o) and absence of transition for v = 4.2 (A) (/ = 30 
Hz) . (b) Autocorrelation of peak pressure versus v illustrating the independence of 
grain mobility on D for D and A^ of 0.66 mm, 9 (o); 0.66 mm, 4 (x); 0.33 mm, 9 
(O); 0.33 mm, 4 (A); 0.46 mm lead, 3 (+); and 0.17 mm steel, 9 (*). 

are independent of D. Quantitatively, the value of Ygm for N < 12 obtained 
from the dispersion data (~ 2.6) is somewhat smaller that the value obtained 
from the local pressure fluctuations (~ 3). A possible explanation for this 
discrepancy is the layer "freezes" from the bottom up as v is reduced. For 
2.6 < V < 3.0 the upper portion of the layer may be mobile but the associated 
mass variations are not seen in the pressure fluctuations because the variation 
in P is small relative to the inherent pressure noise associated with a flat layer 
and/or stress chains in the frozen portion of the layer wash out the pressure 
fluctuations by distributing them over an area comparable to or larger than 
A. 



7 Conclusion 



Two length scales, the layer depth H and the particle diameter D, control 
the scaling properties of granular waves and granular wave patterns. As a 
guide. Table 2 provides a summary of the relevant dimensionless parameters 
and their associated transition values and figures. The horizontal mobility of 
grains appears to be determined by the container velocity — above a critical 
value ~ grains are mobile while below Vgm, the layer depth remains 

uniform throughout the oscillation cycle. The physical mechanism responsible 
for the transition remains to be understood. The sharpness of the grain mobil- 
ity transition stands in contrast to the gradual change in properties of g 
the length scale of disturbances approaches the mean free path. Additionally, 
it is surprising that finite grain size effects are made manifest via v rather 
than H/D or X/D as mentioned in Sec. 5. 
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Table 2 

Summary of parameters, transitions and associated figures. 

Symbol Description Value Figures 



r 


Dimensionless Acceleration 








Flat-layer /Waves Transition 


2.2 (r i) 


3-5 


r 


Dimensionless Frequency 


fVH/g 


5, 8, 9 


/•* 

J ss 


Square/Stripe Transition 


0.33 


1, 6,7 


V 


Dimensionless Velocity 


27rA//v^ 






Grain Mobility Transition 


3{N < 13) 


10-15 



For V > Vgm, the layer behaves as a continuum with A and / scaling with 
H and ^H/g, respectively. With this scaling, dispersion curves for differ- 
ent layer depth and particle size collapse onto the curve given by Eq. 2. 
The square/stripe transition for most parameter values occurs within the 
continuum regime (v > v^^), depends simply on the layer depth H (/*^ = 

fss\jH/g fti 0.33), and shows no systematic dependence on F, A^, D, or parti- 
cle composition (for lead, bronze, and steel). We do not yet understand why 
f*g — 0.33, but speculate it is perhaps related to the dynamic angle of repose. 
For V < Vgm, scaling with H fails. Within this regime, local pressure measure- 
ments support the idea that waves result from a bending of the grain layer 
(see Mujica et al. Ref. [20]), although there is a large discrepancy between 
our value of ~ 3 and that found by Mujica et al. — Vg^ ~ 0.6. Possibly, 
our measurements indicate where horizontal grain mobility stops whereas the 
value reported by Mujica et al. represents the cessation of relative vertical mo- 
tion as well. Another interesting question is the nature of stripe patterns above 
the grain mobility transition. Does wave scaling with H exist in this regime? 
Casual observations show that stripes above the transition have shorter spatial 
and temporal correlations than do stripes in the sloshing regime. 

Finally, the flat-layer/waves transition occurs when the layer remains dilated 
after making contact with the plate. For decreasing F, this transition occurs 
at a nearly constant value of F = 2.2, which is independent of / and N. Asso- 
ciated with the onset of waves is a sudden decrease in the coUision pressure. 
Possibly the least understood aspect of the flat-layer/waves transition is the 
causative connection between the dilation and the pattern. 

There remain numerous unanswered questions concerning wave patterns in 
vertically oscillated granular layers. Increased computational power, better 
simulation techniques, improved measurement capabilities, and application of 
Navier-Stokes-like continuum equations [30] promise significant new insights 
into this system in particular and into the dynamics of granular media in 
general. 
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